Wavefronts may move upstream in semiconductor superlattices 
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In weakly coupled, current biased, doped semiconductor superlattices, domain walls may move 
upstream against the flow of electrons. For appropriate doping values, a domain wall separating 
two electric field domains moves downstream below a first critical current, it remains stationary 
between this value and a second critical current, and it moves upstream above. These conclusions 
are reached by using a comparison principle to analyze a discrete drift-diffusion model, and validated 
by numerical simulations. Possible experimental realizations are suggested. 
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I. INTRODUCTION 



Current instabilities in doped semiconductor superlat- 
tices (SL) have been an active subject of research during 
this decade. For strongly coupled SL, Bloch oscillations 
and Wannier-Stark hopping Q produce negative 
differential conductivity (NDC) at high electric fields. 
This may result in self-sustained oscillations of the cur- 
rent due to recycling of charge dipole domains as in the 
Gunn effect of bulk n-GaAs ||,|| . For weakly coupled SL, 
sequential tunneling is the main mechanism of vertical 
transport. Under dc voltage bias conditions, stationary 
electric field domains may form if doping is large enough 
HQ. Below a critical doping value, the existing charge 
inside the SL may not be able to pin domain walls, and 
current self-oscillations appear [||||. These oscillations 
may be due to recycling of charge monopoles (domain 
walls) or dipoles depending on the boundary condition 
at the injecting contact region (in a typical n + -n-n + con- 
figuration with the SL imbedded between highly doped 
regions, the doping at the emitter region is crucial) pOl . 
Driven chaotic oscillations have also been predicted |j I 
and observed in experiments fl^| . Lastly, there are ways 
to tune the charge inside the SL (and therefore obtain 
stationary domains or self-oscillations) without replacing 
it by a different one. For example, by applying a trans- 
verse magnetic field or by photoexciting the SL jlJJ . 

Transport in weakly coupled SL can be described by 
simple rate equation models for electron densities and 
average fields in the wells, |l5|-|l8|. Many of the effects 
related above have been explained by means of a simple 
discrete drift model jl6 17 l|j2^] . In this model, the tun- 
neling current between two adjacent wells, Ji—n+x, equals 
the 2D electron charge density at well i times a drift ve- 
locity, which depends on the electric field at the same 
well. By starting from a microscopic sequential tunneling 



model, it has been shown that the discrete drift model is 
a good approximation at low temperatures and for fields 
above the first plateau of the SL current- voltage charac- 
teristic Jl^,|2l[ . For low dc voltages on the first plateau, 
a discrete diffusion (which is a nonlinear function of the 
field) should be added. This term contains the contri- 
bution to Ji^i+i of the tunneling from well i + 1 back 
to well i (which vanishes for large enough electric fields) 
p8| , pl| . In this paper we report an interesting conse- 
quence of electron diffusivity at low fields: if the current 
is sufficiently high, and so is the doping, a domain wall 
(monopole wave) which connects two domains may travel 
in a direction opposite to the flow direction for electrons 
(i.e., upstream, in the positive current direction!). This 
striking phenomenon is contrary to the usual situation: 
a monopole either moves downstream (in the direction of 
the flow of electrons), or it remains stationary, jt9). We 
substantiate our claim both by numerical simulations of 
the discrete drift-diffusion model and by rigorous mathe- 
matical analysis based upon a comparison principle p2[ . 
Mathematical analysis yields useful bounds for critical 
values of current and well doping, and for monopole ve- 
locity. 

There are related fields for which differential-difference 
equations (similar to discrete drift-diffusion models) 
model the systems of interest. Well-known are propa- 
gation of nerve impulses along myelinated fibers, mod- 
elled by discrete FitzHugh-Nagumo equations p3| , p4| , 
motion of dislocations (2^j2(| and sliding charge density 
waves [ 27| , modelled by variants of the Frenkel-Kontorova 
model Jp8| , etc. The theory of wavefront propagation has 
been developed for some of these models, which are sim- 
pler than ours: convection is typically absent from them 
and diffusion is purely linear 

The rest of the paper is as follows. We write the drift- 
diffusion model with appropriate boundary conditions in 
Section ||. There we render these equations dimension- 
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less and explain the results of numerical simulations on 
a current biased infinitely long SL. Furthermore, we find 
by numerical simulations that our results for infinite SL 
may be realized in finite SL with appropriate boundary 
conditions under constant current bias. The theoretical 
analysis based on the comparison principle is presented in 
Section IV contains our conclusions. Finally 



,(/) 



Section [II 



some material of a more technical nature is relegated to 
the Appendices. 



II. DISCRETE DRIFT-DIFFUSION MODEL 

A. Equations and boundary conditions 

At low enough temperatures (much less than a typi- 
cal Fermi energy of a SL well measured from the first 
subband, say 20 meV or 232 K), the following discrete 
drift-diffusion equations model sequential vertical trans- 
port in a weakly doped SL Jf^,^]} : 

e d Fi , ni v(Fi) D{ n l+1 -n ; = 



e dt d + w 



(d + wf 

F i -F i _ 1 = ^(n i -N%), (2) 



Eq. (|lj) is Ampere's law establishig that the total current 
density, eJ, is sum of displacement and tunneling cur- 
rents. The latter consists of a drift term, en,iv(Fi)/(d + 
w), and a diffusion term, eD(Fi) (rij+i — rii)/(d + w) 2 . 
We have adopted the convention (usual in this field) that 
the current density has the same direction as the flow of 
electrons. Eq. (g) holds for i = 1, . . . , JV - 1. Eq. (§) is 
the Poisson equation, and it holds for i = 1, . . . , N. n,; is 
the 2D electron number density at well i, which is singu- 
larly concentrated on a plane located at the end of the 
well. Fi is minus an average electric field on a SL period 
comprising the ith well and the ith barrier (well i lies 
between barriers i — 1 and i: barriers and N separate 
the SL from the emitter and collector contact regions, 
respectively). Parameters e, d, w, and Np are well per- 
mittivity, barrier width, well width and 2D doping in the 
wells, respectively. 

Drift velocity and diffusion coefficient are depicted in 
Fig. for the 9nmGaAs/4nmAlAs SL of Ref. @. We have 
obtained them from microscopic calculations presented in 
Rcf. (which is appropriate for these sample param- 
eters^) by setting v(F) = J(N%,N%,F) (d + w)/N l £ 
and D(F) = -{dJ{N l £, N%, F)/dn l+1 )(d + w) 2 . Here 
e J (rii, n.j+1, Fi) is the tunneling current between wells i 
and i + 1, Jj^+i. We assume that the tunneling current 
is a function of the average field at the ith SL period, 
Fi = F, and of the 2D electron densities at wells i and 
i + 1, rii and nj+i, respectively. Notice that our model 
for the tunneling current, 



iJ{m,ni+\,Fi) = 



en i v(F l ) 
d + w 



eD(Fi) 



H+l 



(Fi) 



+iv^(F i ) 



w 



(3) 



is reasonable for temperatures much lower than a typ- 
ical Fermi energy in the wells measured from the first 
subband (say 20 meV), p3j. The tunneling current den- 
sity should change sign if we reverse the electric field 
and exchange the electron densities at wells i and i + 
J(rij, rij+i, Fi) — — J(m+i, rii, —Fi). This inversion sym- 
metry implies 

v {f) (-F) = v {b) (F) and v(-F) = -v(F), 

where v^(F) = D(F)/{d + w) and v^(F) = v(F) + 



(F). See Figure Md) 



jib) 

Equations (|lj) and"~(|2|) should be supplemented with ap- 
propriate bias, initial and boundary conditions. Among 
possible bias conditions, we shall consider the extreme 
cases of current bias (J(t) specified) and voltage bias: 



N 



(d + w)Y,F i = V. 



(4) 



with specified V = V(t). Using (f§) ignores potential 
drops at the contact regions and at barrier 0, and it 
overestimates the contribution of barrier N by a factor 
1 + w/d pl[ |. These contributions are negligible for long 
SL (N — 40 or larger), so that we shall adopt the simpler 
expression (||). Appropriate boundary conditions have 
been derived under the same approximations as in (|l|) 
|l|. They are 



e dt d + w 



e dF N n N wW(F N ) 



dt 



d + w 



J(t), 



(5) 
(6) 



where the emitter current density, eje (F), the emitter 
backward velocity, w^'(F), and the collector forward ve- 
locity, w^\F) are functions of the electric field depicted 
in Fig. 3 of Ref. |2l]] for contact regions similar to those 
used in experiments § . 

To analyze the discrete drift-diffusion model, it is 
convenient to render all equations dimensionless. Let 
v(F) reach its first positive maximum at (Fm,vm)- We 
adopt F M , Np, v M , v M (d + w), eN%VM/(d + w) and 
£F M (d+w)/(eN%v M ) as the units of Fi, n u v(F), D{F), 
eJ and t, respectively. For the first plateau of the 9/4 SL 
of Ref. (|, we find F M = 6.92 kV/cm, N% = 1.5 x I0 11 
cm~ 2 , vm = 156 cm/s, vm (d+w) = 2.03 x 10~ 4 cm 2 /s 
and eN%v M /(d + w) = 2.88 A/cm 2 . The units of current 
and time are 0.326 mA and 2.76 ns, respectively. Then 
(0) to (|) become 



dEi 

—jjr + v ( E i) n i ~ D(Ei) (n i+1 - rii) = J, 



E, 



{d + wf 



Ei-i =v(m- 
1 N 

N ^ 



1). 



(7) 
(8) 

(9) 
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Here we have used the same symbols for dimensional 
and dimensionless quantities except for the electric field 
(F dimensional, E dimensionless). The parameters v — 
eNp/(e F M ) and <j> = V/[F M N(d+w)] are dimensionless 
doping and average electric field (bias), respectively. For 
the 9/4 SL, v » 3. We recall that i = 1, . . . , N - 1 in 
(fjj) and i=l,...,JVin(§). The boundary conditions (j^) 
and M) become 



dE 

—— + J e {E Q ) - w e (E ) m = J, 
at 



dE 



N 



dt 



+ w c (E N )n N = J, 



(10) 
(11) 



where 



Je{E ) 



ji f \F M E ){d + w) 



w e (E ) 
w c (E N ) -- 



N%v M 
w^(F m Eq) 

VM 

w 

VM 



(12) 



Figure || shows J e , w e and w c as functions of the elec- 
tric field. They are dimensionless versions of the curves 
plotted in Figure 3 of Ref. flj. 



B. Numerical simulations 

Simple solutions of the drift-diffusion equations ((?]) - 
(|J) under constant current bias are stationary or mov- 
ing monopole wavefronts connecting two electric field do- 
mains. Let us consider monopole solutions with profiles 
{Ei} which are increasing functions of i, for they are com- 
patible with realistic boundary conditions in which the 
emitter region is highly doped We have simulated 
numerically on a large SL, 



dE, _ D(Ei)+v(Ei 
dt v 
D(Ek) 



(Ei-i — Ei) 



(E i+1 -Ei) = J-v(Ei), 



(13) 



with fixed J, which is equivalent to (j^) - (S). Let 
EW(J) < E@>(J) < EW(J), be the three solutions of 
v(E) = J for v m < J < 1, where (E mi v m ) is the min- 
imum of v{E) for E > 1. For the 9/4 SL of Figjjl, 
E m = 9.8571, v m = 0.02192. We have simulated |yf) 
for different values of v > and of J £ (v m ,l). The 
initial condition was chosen so that Ei — » E^ x \j) as 
i — > -co, and Ei -> E^(J) as i — y oo. We observed 
that, after a short transient, a variety of initial condi- 
tions sharing these features evolved towards either a sta- 
tionary or moving monopole. For systematic numerical 
studies, we therefore adopted an initial step like profile, 
with E % = EW(J) for i < 0, Ei = E^(J) for % > 



and Eq — E^(J). The boundary data were taken to be 
E- N = EW(J), E n = E^(J) with N large. 

Our results show that the dimensionless doping v de- 
termines the type of solution of (13) which is stable. 
There are two important values of v\ <vi- 

• For < v < v\ and each fixed J £ (v m , 1), only 
traveling monopole fronts moving downstream (to 
the right) were observed. For v > v\, stationary 
monopoles were found. According to the arguments 
of Wacker et al for the discrete drift model 
with D(E) = 0, stationary monopoles exist for di- 
mensionless doping larger than a critical value. An 
upper bound for this critical doping is 



En 



1 



(14) 



which equals v c = 0.198 for our numerical example. 
We have found that V\ = 0.16. This agreement 
with results obtained assuming D(E) = is not 
surprising: we shall prove in Section |l| that ( |l4| ) 
holds as well for the model ([?]) - ^ with nonzero 
diffusivity. 

• For v\ < v < i/2, traveling fronts moving down- 
stream exist only if J £ (v m , Ji(i/)), where Ji(v) < 
1 is a critical value of the current. If J £ ( Ji(i'), 1), 
the stable solutions are steady fronts (stationary 
monopoles). We have found that v 2 = 0.33. 

• New solutions are observed for v > As be- 
fore, there are traveling fronts moving downstream 
if J £ {v m , J\(y)), and stationary monopoles if 
J £ (Ji(u), J2(v)), Jziv) < 1 is a new critical cur- 
rent. For J<z{v) < J < 1, the stable solutions of ( |l3| ) 
are monopoles traveling upstream (to the left). As 
v increases, J\(y) and Ji(y) approach v m and 1, 
respectively. Thus stationary solutions are found 
for most values of J if v is large enough. 

Figure || depicts J\(v) and Jiiy) as functions of v. 
Notice that J% decreases from J\ = 1 to J\ = v m as v 
increases from v\. Similarly, J 2 decreases from J 2 = 1 to 
a minimum value J 2 w 0.53 and then increases back to 
J 2 — 1 as v increases. Monopole velocity as a function of 
current has been depicted in Figure ^ for four different 
doping values, v = 0.5, v = 1, v = 3 and v = 10. For 
larger v, the interval of J for which stationary solutions 
exist becomes wider again, trying to span the whole in- 
terval (v m , 1) as v — * oo. For very large v, the velocities 
of downstream and upstream moving monopoles become 
extremely small in absolute value. 

Notice that if we use the complete sequential tunnel- 
ing current instead of the drift-diffusion approximation 
(||) in Eq. |l]), the situation is the same. Figure || de- 
picts monopole velocity versus current for well doping 
corresponding to the 9/4 SL of Ref. ||. Results obtained 
with the complete sequential tunneling current or with 
approximation (||) (corresponding to Fig. [| with v = 3) 
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are compared. Both velocity curves are similar, and their 
quantitative discrepancies are irrelevant in view of the 
uncertainties involved in a theoretical calculation of the 
tunneling current (typically the off-resonance current is 
larger than the theoretical prediction). 

Once different stable monopole solutions (moving ei- 
ther downstream or upstream, stationary) have been 
identified, we raise the natural question of whether they 
are compatible with boundary conditions. Another series 
of numerical simulations was carried out to answer this. 
We solved numerically ( |l3|) for a current biased finite SL 
(N = 40) with boundary conditions (jlO) - Our re- 

sults are depicted in Figure || for realistic doping at the 
contact layers. We observe that the emitter boundary 
condition results in the creation of a charge accumula- 
tion layer near this contact. A charge depletion layer is 
formed near the collector contact as a result of the cor- 
responding boundary condition. Except for these layers, 
existence and configuration of monopoles moving down- 
stream, upstream or remaining stationary, agrees with 
the previous simulations (corresponding to an infinitely 
long current-biased SL with a monopole-like initial con- 
dition) . 



III. MATHEMATICAL ANALYSIS OF 
TRAVELING MONOPOLES AND STATIONARY 
SOLUTIONS 



In this Section, we study theoretically moving or sta- 
tionary monopoles on an infinitely long, current-biased 
SL. Our findings will confirm the picture suggested by 
the numerical simulations of the previous Section for 
any doped weakly coupled SL. Furthermore, we shall 
prove stability of the different monopole solutions and 
find bounds for the critical values of v and Ji. Our 
results are based upon and extend ideas first proposed 
by J. P. Keener for discrete FitzHugh-Nagumo equations, 
corresponding to signal transmision in myelinated neu- 
rons 124] . Mathematically analogous problems arise in 
models of propagation of defects in crystals (2f| . These 
problems have the following structure, 

dE- 

~^-d (E l+1 - 2E { + = J - v(Ei), (15) 

which is much simpler than (|l3|). Here the parameter 
d > is a constant diffusion coefficient, and v{E) a 'cu- 
bic' function with three branches as the electron drift 
velocity of Fig. [TJ. 

For (15), there are critical values of J, J\ and J2, 
characterizing wavefront behavior |p4[ . For J > J 2 (d) , 
there exist wavefront solutions of Jig) moving upstream 
(to the left). For J < Ji(d), there are wavefronts mov- 
ing downstream (to the right), whereas for J\(d) < J < 
J2{d), stationary fronts exist. The width of the interval 
(Ji(d), Jiid)) is an increasing function of d. 



A. Propagation failure and stationary solutions 

In Appendix ^ we state and prove a comparison prin- 
ciple for (§J). As a consequence, if our initial field profile 
is monopole-like [monotone increasing with well index, 
and sandwiched between E^(J) and E^ 3 \j)], so is the 
electric field profile for any later time t > 0; see Appendix 

{Ei(0)} increasing with i ===> 
E^(J) < Ei{t) < E l+1 {t) < E^(J), Vi, t > 0. 

We now obtain sufficient conditions for an initial 
monopole not to propagate upstream or downstream. 
Under these conditions, the monopole may remain sta- 
tionary or move downstream or upstream, respectively. 
Let us start with a condition pinning the left tail of a 
monopole. As i^-i < £7, and E i+ i < E^ 3 '(J), we have: 



dEj 
dt 



D{Ei 



(Ei+i — Ei) 



D(Ei)+v(Ei) 



(Ei-i - Ei) + J — v(Ei) 
-Ei]+J- v(Ei) < 0, 



provided there exist a; < 6; such that 

W-[E-EW(J)\ > J-v(E), Ee{ai,b{), (16) 
v 

and then we choose some initial field, Ei(Q) £ (a/, bi). 
The previous inequality then implies E^t) G (E^ (J), b t ) 
for all t > 0. This in turn forbids a monopole to move 
upstream (to the left). We say that condition ( [La ) pins 
the left tail of the monopole. Whether such (a/, bi) exist, 
depends on the parameters v and J; see Figure [?]. 

Let us now pin the right tail of a monopole. As 
E i+1 > E l and E l ^ 1 > E^(J), we have 



dE, 
~dt 



v 

D{Ei) 



v{Ei) 



> 



D(E. t ) + v(E l 



{E^ - Ei) + J - v(Ei) 
{E^(J)-Ei] + J-v(Ei)>0, 



provided there exist a r < b r such that 

[E -E*\j)]<J-v { E), 

V 

E e (a r ,b r ), 



(17) 



and we choose some initial field, Ei(0) 6 (a r ,b r ). The 
previous inequality then implies Ei(t) G (a r , E^ ( J)) for 
alii > 0. A monopole cannot then move downstream (to 
the right), and we say that its right tail is pinned. Figure 
@ illustrates our arguments: for fields larger than a r , the 



4 



Ei's tend to increase above a r toward E^(J). Then the 
monopole cannot move downstream. For Ei < b\, the 
fields tend to E^\j), and the monopole cannot move 
upstream. As before, the existence of (a r , b r ) depends on 
the values of v and J; see Figure 0. 

Figures @ show the curves J-v(E), D(E) (E-E^)/u 
and [D(E) + v(E)} (E - E^)/u for v = 3 and differ- 
ent values of J. At J = 0.08, Figure 0(a) shows that 
there is an interval (ai, bi) as in (|l6|), but no interval 
(a r , b r ) as in jl7| ) exist. Then the left tail of a monopole 
is pinned, but its right tail is free. In this conditions, 
a monopole may move downstream. Figure 0(b) shows 
a monopole with both its left and right tails pinned for 
J = 0.2. Then our theory implies that wavefront prop- 
agation fails and a monopole-like stationary solution is 
stable. Numerical simulations show that there are sta- 
tionary solutions when J G [0.09,0.53]. Finally, Figure 
0(c) shows that, if J = 0.34, the right tail of a monopole 
is pinned, but not its left tail. Under these conditions 
a monopole may move upstream. For larger v, the es- 
timates become sharper. For instance, when v = 10, 
@ and @ hold for J G [0.05,0.45]. Direct numeri- 
cal simulations show that stationary solutions exist for 
J G [0.04,0.55]. Systematic use of these criteria allows 
us to estimate the critical doping values Vj and critical 
current values J^v)^ = 1,2 defined in the previous Sec- 
tion; see Appendix [b|. Instead of looking for J\ iy) and 
Ji(y), it is more convenient to look for their inverse func- 
tions, which we may call v\(J) and v^(J). According to 
Figure |, the inverse function of J2{ v ) is two- valued, and 
its two branches are z^(J) < v i{J)- We have found the 
following upper bounds vib{J) and ^(J) for v\{J) and 
i>2 (J), respectively: 



E m -EW(J) 



J — V 



m, 



l - J 



(18) 
(19) 



If v > V\b(J), the right tail of the monopole is pinned, 
whereas the left tail of the monopole is pinned if v > 
V2b(J)> see Appendix [^. Notice that V\b{J) is a de- 
creasing function of J. Therefore the critical value V\ 
(above which there are stationary solutions) is smaller 
than ^ib(l), which is exactly Wacker et al's bound, (14). 
This explains why the bound (|l4|) gives surprisingly good 
results even for the first plateau of the SL current-voltage 
characteristics [despite having been obtained under the 
assumption D(E) = 0] |ll],[l8|]. Notice that the bound 
( |l9| ) is reasonable for large dopings and currents J « 1. 
See Figure^ for a comparison between the critical curves 
J\(v) and J%(v) and the bounds (Il8j) and (|l9|). 



B. Propagation: traveling fronts 

Having shown that only one tail of a monopole is 
pinned suggests that the monopole may move in the op- 



posite direction. Direct simulations show that this is of- 
ten the case, and we will prove this now. 

An upstream traveling wave solution of (|l3| ) may have 
the form 



Ei(t) = w(i + ct), c> 0. 



(20) 



We will look for an electric field profile w(z), z = i + ct, 
which is not an exact solution of ([l^) , but instead it sat- 
isfies 

dw D(w(z)) + v(w(z)) 
c — [w(z — 1) — w{z)\ 



D(w(z)) 



■){z + 1) - w(z)} + v(w(z)) - J < 0. (21) 



If this subsolution is initially below an initial field profile, 
i.e. w(i) < Ei(G), for all i, then the comparison theo- 
rem of Appendix [A] guarantees that Ei(t) > w(i + ct) for 
later times. As w(i + ct) moves upstream, so does Ei(t), 
and the electric field profile corresponds to a monopole 
moving upstream with velocity at least c. See Figure 0: 
a subsolution "pushes" the monopole upstream, whereas 
a supersolution (defined below) "pushes" the monopole 
downstream. 

How do we find a reasonable subsolution? An idea 
is to try a piecewise continuous solution which equals 
E^'(J) for z < zq and a larger constant A, A G 
(E^(J),E^(J)) [and therefore v(A) — J < 0], for 
z > z\, with z\ > zq. For zq < z < z±, w{z) is an unspec- 
ified smooth increasing function with w{zq) — E^( J) 
and w(zi) = A. Now we shall select conveniently the 
numbers zq, z\, c and A so that ( |2l| ) holds. Clearly, ( |2l| ) 
holds for z + 1 < zq and for z — 1 > Z\. Suppose that 
< Z\ — zq < 1. Then there are five possibilities: 

1. z < zq, z + 1 > z\. Then w(z — 1) = w{z) = 
E {1 \j), 10(21+1) = A, which inserted in @ yields 
-D(E^) (A - E<V)/v < (obviously true). 

2. z - 1 < z Q , z > zi. Then w(z - 1) = E^\j), 
w(z) — w(z + 1) = A, which inserted in (H|) yields 



J - v(A) > 



D(A)+v(A). 



[A-E (1 \J)}. (22) 



3. Zq < z - 1 < Zf Then E^{J) < w(z - 1) < A 
and w{z) = w(z + 1) = A, which yields J — v(A) > 
[D(A) +vCA)] [A - w(z - l)]/v. This inequality 
holds if (|22| ) does. 

4. zq < z < z x . Then w(z - 1) = E^{J), BW(J)< 
w{z) < A and w(z + 1) = A. Inserting this in (|2l|), 
we find 

c < D{w{z)) + v{v}{z)) (1) _ 
dz 



v 



D{w{z)) 



[A - w(z)] + J - v(w(z)). 
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Let us now assume that we can select A G 
(£(2)^(3)) suc h that the right hand side of this 
expression is positive, say 

D(w)+v(w) ^(i) (J) _ w] + D H (A-w) 

V V 

+ J - v(w) > 6 > 0, E {1) < w < A, (23) 

and that we choose c so that cdw/dz < S. Then 
© holds. 

5. z Q < z + 1 < z\. Then w(z - 1) = w(z) = E^ (J), 
and E^{J) < w(z + l) < A, which inserted in ( |2l| ) 
yields -D(E^) [w(z + 1) - E^]/v < (obviously 
true). 



Summarizing the previous arguments, provided ( |22| ) and 
( p3| ) hold, w(z) is a subsolution obeying (pl|). The param- 
eter A can be found graphically. First of all, we depict 
the functions J-v(E) and f x (E; J) = [D{E) + v(E)] [E- 
E^\j)]/v. Possible values of A are those E for which 
J — v(E) > /i(-E; J). For such A, we may plot the left 
side of (p3j), 



f2(E;J,A) 



_ D(E) + v(E) 



[EU(J)-E] 



D(E)(A-E) 
+ — - + J- v(E). (24) 



Ef 2 (E;J,A) >0for£e (E^(J),A), then the selected 
value of A allows us to construct the sought subsolution. 
See Figure |l0| for a practical realization of this graphical 
construction. 

We have proved rigorously that monopoles may 
move upstream under favorable circumstances. Our 
proof using subsolutions may yield a very practical 
additional bonus: an upper bound, c* , for the ve- 
locity of the monopole. Let us choose 5( J, A) = 
™n E( i )<E<A f 2 (E;J,A), z x - z = 1, and w(z) = 
[A - E^(J)](z - z ) for z < z < Zl . Then c* = 
6(J,A)/[A - EW(J)\. In Figure |, — c* is represented 
by a line of thick dots for doping v = 3 corresponding to 
the 9/4 SL. 

In a similar vein, we can construct supers olutions which 
push the monopole field profile to the right; see Figure 
|9j. Now we start from a monopole profile moving down- 
stream, 

Ei(t) = w(i-ct), c> 0. (25) 
The electric field profile w(z), z = i — ct should satisfy 
dw D(w(z)) + v(w(z)) 



dz v 
D{w{z)) 



[w(z — 1) — w(z) 



[w(z + l)-w{z)] + J-v(w(z)) <0. (26) 

We seek a piecewise continuous supersolution which 
equals a constant, A, A £ (E^ (J) , (J)) , for z < z , 



and w(z) — E^(J) for z > z\, with z\ > z . For 
Zq < z < Zi, w(z) is an unspecified smooth increasing 
function with w(zq) — A, and w(zi) — E^(J). As 
for subsolutions, we now select conveniently the num- 
bers zq, zi, c and A so that (26) holds. Clearly, ( |2^ ) 
holds for z + I < zq and for z — 1 > Z\. Suppose that 
< Z\ — zo < 1. An analysis of the remaining five possi- 
bilities yields the following criteria to hold for w(i — ct) 
to be a supersolution: 



j _ v{A) < _£!£. [S (3) (J) _ A] = h(A . J)) 

/ 4 (w; J, A) = (w - A) 

[£ (3) (J) - w] + J - «(«>) < -5, 



+ 



for 



A < w < E {3) (J), 

C — < 0. 



(27) 

(28) 
(29) 



Provided such w(i — ct) is found, solutions Ei(t) of 
with Ei(0) < w(i + t) will satisfy E t (t) < w(i — ct + r) 
and propagate to the right with speed larger than c. r is 
a constant which can be conveniently chosen to keep the 
monopole profile below the supersolution. Figure [ll] il- 
lustrates the graphical construction of the supersolution 
by checking that ( p7j ) and ( p8| ) hold for particular values 
of J and v. 

As in the subsolution case, an upper bound c* for the 
monopole velocity c is estimated by choosing —S( J, A) — 
max A<E<E (3)f 4 (E;J,A), z\ - z = 1, and w(z) = 
[EW(J) - A](z - z ) for z < z < z x . Then c* = 
5 (J, A)/[E( 3 '(J) - A]. In Figure |, c* is represented by a 
line of thick dots for doping v — 3 corresponding to the 
9/4 SL. 

We can now summarize the results obtained from sub 
and supersolutions; see Figures and 0. We find rea- 
sonably good upper bounds for the absolute value of 
monopole velocity. Furthermore, for v — 3, conditions 
(^||) and ( p3|) hold for J = 0.6 and A = 12, whereas 
conditions ^2]) and (§|) hold for J = 0.05 and A = 3. 
Therefore, monopoles move downstream for J < 0.05 and 
they move upstream for J > 0.6. Direct numerical simu- 
lations show that: (i) the estimate J\ = 0.05 for the first 
critical current can be improved to J\ = 0.08; and (ii) 
J2 = 0.6 for the second critical current can be improved 
to J 2 = 0.54. 



IV. CONCLUSIONS AND FINAL COMMENTS 

We have presented a theory of monopoles moving 
downstream or upstream on an infinitely long doped, cur- 
rent biased superlattice when the fields are on the first 
plateau of the current -voltage characteristic. This theory 
has been corroborated with numerical evidence, which 
sharpens our results. Furthermore, we have simulated 
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a 40- well 9nmGaAs/4nmAlAs SL M under do ping and 
contact conditions similar to experimental ones 21 1, but 
under constant current bias conditions. This situation is 
different from the usual case of voltage bias conditions. 
We have obtained that it is possible to observe monopole 
wavefronts moving upstream when the current is kept at 
large enough levels. Together with our theoretical bounds 
for critical currents and dopings, this numerical predic- 
tion could be used to set up an experiment to observe 
this striking phenomenon. For this purpose, we would 
need an initial condition corresponding to a monopole 
separating two electric field domains at high enough cur- 
rent. In an ideal world, this situation could be obtained 
by first fixing a low dc voltage for the 9/4 sample at 
a value near the top of one of the first branches of the 
current - voltage characteristics. Then we could switch 
from voltage to current bias conditions. The outcome 
would be a monopole moving upstream until the emitter 
region is reached. Presumably an idea of the field distri- 
bution corresponding to this situation could be obtained 
by time-resolved photoluminescence measurements |l. 

There are technical problems that must be overcome if 
one wants to observe these features in real experiments: 
when we switch, there will always be a Faraday-like in- 
ductive pulse which will probably perturb the state of the 
system in an uncontrolled way. There are other possible 
biases we could think of. Under dc voltage bias, upstream 
moving monopoles are probably created for a short time 
during relocation experiment pc| . In these experiments, 
one has a doped SL with a current-voltage characteris- 
tics correponding to multiple stationary monopole solu- 
tion branches. Voltage is set at a particular value near 
the end of a branch, so that the field profile is that of a 
monopole layer connecting a low to a high field domain. 
Let the monopole layer be located at well i (counted 
from the emitter contact). Then the voltage is suddenly 
and appropriately increased. After a certain time, the 
field profile settles to a new situation corresponding to 
a monopole layer centered at well i — 1 J3(J . This could 
be an indication of a monopole moving upstream, albeit 
for a short time. To increase this time, we could try to 
set a hybrid bias (between current and voltage bias) by 
including a finite series resistance in our external circuit. 
Additional theoretical and numerical work is needed to 
explore these possibilities. 
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APPENDIX A: COMPARISON PRINCIPLE 

The ma in t heorem which we use to prove our results 
in Section III is the following comparison principle: 
Theorem A.l Let Ui(t) and Li(t), i G Z, be differen- 
tiable sequences such that 

^-d 1 (U i )[U i+1 -U i ] 

-thWpi-i-Ud- f(Ui)> 

— — d\ (Li) [Li+i — Li] 

-d 2 (Li)[Li-i - Li] - f(Li), 
Ui(0) > Li(0). 



(Al) 



where f, d\ > and d 2 > are Lipschitz continuous 
functions. Then, 

Ui(t)>Li(t), t>0,ieZ 



In our discrete drift-diffusion model, 

V V 

f(E) = J-v(E). 

Proof: The proof is by contradiction. Set Wi(t) = 
Ui(t) - Li(t). At t = 0, Wi(0) > for all i. Let us 
assume that Wi changes sign after a certain minimum 
time ti > 0, at some value of i, i = k. Thus Wk(t\) = 
and dWk/dt < 0, as t — ► t\. We shall show that this 
is contradictory. At t = ii, there must be an index m 
(equal or different from k) such that W m (ti) = 0, while 
its next neighbor W m +j(ti) > (j is either 1 or -1), and 
Wi(t±) = for all indices between k and m. For other- 



wise Wk should be identically for all k. Equation (Al) 
implies 



dW„ 
dt 



■(h) > di(am(*i))wWi(*i) 

+d 2 (U m (h))W m - 1 (t 1 ) > 0. 



This contradicts the fact that dW m /dt should have been 
nonpositive as t — » ti, for W m (ti) to have become zero 
in the first place. 



Corollary A.l Any solution Ei(t) of ( ]i4 ) with ini- 
tial data Ei(Q) S (£ (1) (J), E^(J)) satisfies E 2 (t) G 
(EW(J),EW(J)) fort>0. 

Proof: Apply Theorem A.l first with L, = E^(J) 
and Ui = Ei, then with L; = E % and U t = E^(J). 

Corollary A. 2 If Ei(0) is monotone increasing, that is, 
Ei(0) < Ei+\(0), then, Ei(t) is also monotone increas- 
ing, i.e., Ei(t) < Ei + i(t) for t > 0. 
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Proof: Apply Theorem A.l with Li = Ei(t) and 
Ui=E i+1 (t). 

Remark. Strict inequalities in these theorems can be re- 
placed by inequalities and the corresponding statements 
still hold. However the proofs become rather more tech- 
nical and involved. 



APPENDIX B: BOUNDS FOR CRITICAL 
DOPING VALUES 

We want to estimate the curves v\b{J) and v^ b {J) de- 
To estimate ia^(J), assume that 
and v is large. The left tail of a monopole is 



fined in Section [III 
J — ► 1- 

pinned if Eq. (|n]) holds. For large currents, ([L6|) cer- 
tainly holds if the curve corresponding to the left side of 
the inequality is above that of the right hand side, for 
E = 1 (this is possible because D(E) decreases rapidly 
to zero as the field increases). Setting E = 1 in (|16|), we 
obtain 



£(!)[! -£( 3 ) (J)] 



> J - 1. 



In turn, this implies v > ^i(J), defined in (fl9l). This ar- 
gument fails for the small values of J used to draw Figure 
0. We believe that quite different reasoning is needed to 
estimate z/^(J). 



The same argument yields our estimate vu(J) of (18). 
For Jl7| ) to hold, the curve corresponding to the left side 
of the inequality should be below that of the right hand 
side for E — E m . As D(E m ) « 0, we obtain 



v(E m ) 



[E m -E^{J)\ < J-v(E m ), 



which yields (jig). Fig. || shows that the bound (18) is 
reasonably good for all eligible values of v and J. 
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FIG. 1. (a) Drift velocity and (b) diffusion coefficient 
for a 9nmGaAs/4nmAlAs SL. Doping at the wells is 



cm , whereas at the contact regions, 

-3 



Nd — 2 x 10 cm . (c) Dimensionless drift velocity, v(E), 
diffusion (equivalent to backward tunneling velocity), D(E), 
and forward tunneling velocity, v^\E) = v(E) + D(E). (d) 
Extension of the dimensionless diffusivity to negative values 
of field. We have D{-E) = « (/) (£:). The same formula yields 
the extension of v^'(E) to negative fields. Then v(E) is an 
odd function of E. 
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FIG. 3. Critical currents Ji and J2 as functions of the 
dimensionless doping v. Monopoles move downstream for 
v m < J < J\(v), are stationary for Ji(f) < J < Ji(y), and 
move upstream for Ji(y) < J < 1. Dashed lines in this figure 
represent the bounds ^ib(J) and u^ b {J). 
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FIG. 2. Dimensionless functions of the electric field for the 
contact regions, (a) Current at the emitter, J e (E). (b) Back- 
ward velocity at the emitter, w e (E). (c) Forward velocity at 
the collector, w c (E). 



FIG. 4. Velocity of a monopole wavefront as a function of 
J for four doping values v — 0.5, 1, 3, and 10. Monopoles 
with negative velocity move upstream. For doping v = 3 cor- 
responding to the 9/4 SL, we have also represented bounds 
for the velocity as lines of thick dots. 
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FIG. 5. Monopole velocity as a function of current for the 
9/4 SL. Comparison of results for the discrete drift-diffusion 
model (1) - (2) with v = 3 (thick line), and those obtained by 
using the exact tunneling current eJ(n,i, ni+i, Fi) (thin line) 
instead of approximation (3) . 
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FIG. 6. Numerical simulations of the drift-diffusion model 
with realistic boundary conditions at the contact regions, (a) 
Monopole moving downstream for J = 0.023, (b) stationary 
monopole for J = 0.3, and (c) monopole moving upstream for 
J = 0.9. In all cases, diamonds correspond to the profile at 
t = and squares to the profile at the largest positive time. 




11 















0.6 
0.4 


i 


i 


i i 


1 1 ; 
(C)j 


i 




0.2 

















- L- 






- — 4 




















-0.2 














-0.4 














-0.6 






i i 


i i 


i 




2 


4 


6 8 
E 


10 12 


14 



FIG. 7. (a) Functions of the electric field establishing pin- 
ning of the monopole tails for v = 3, J = 0.08. Solid line: 
J - v(E), dashed line: D(E) [E - E (3 \j)]/v, dotted line: 
[D(E) + v{E)] [E - £«(J)]/ia (b) Same plots as in (a) for 
v — 3, J = 0.2. (c) Same plots as in (a) for v = 3, J = 0.34. 
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FIG. 8. A monopole field profile. The field values bi and 
a r are indicated by horizontal dashed lines. Vertical arrows 
indicate that the field at those SL periods either (i) decrease 
toward E^ (J), and therefore the monopole left tail is pinned; 
or (ii) increase toward E^(J), and therefore the monopole 
right tail is pinned. 



FIG. 9. A monopole field profile, sub and supersolutions. 
The supersolution is always above the real values of the field, 
and therefore it pushes the monopole to the right. The sub- 
solution pushes (from below) the monopole to the left. 




FIG. 10. Curves determining the subsolution for v = 3, 
J = 0.6, A = 12. Solid line: J - v(E); dashed line: J); 
dotted line: f2(E; J, A). 





FIG. 11. Curves determining the supersolution for v — 3, 
J = 0.05, A = 3. Solid line: J - v(E); dashed line: f 3 (E; J); 
dotted line: f 4 {E; J, A). 
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